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Gauge-invariant observables for quantum gravity are described, with explicit constructions given 
primarily to leading order in Newton’s constant, analogous to and extending constructions first given 
by Dirac in quantum electrodynamics. These can be thought of as operators that create a particle, 
together with its inseparable gravitational field, and reduce to usual field operators of quantum 
field theory in the weak-gravity limit; they include both Wilson-line operators, and those creating a 
Coulombic field configuration. We also describe operators creating the field of a particle in motion; 
as in the electromagnetic case, these are expected to help address infrared problems. An important 
characteristic of the quantum theory of gravity is the algebra of its observables. We show that 
the commutators of the simple observables of this paper are nonlocal, with nonlocality becoming 
significant in strong field regions, as predicted previously on general grounds. 
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I. INTRODUCTION 

Gravity is widely believed to be described by a quantum-mechanical theory, rather than one that requires an 
extension or modification of quantum mechanics. If this is the case, the structure of quantum mechanics (suitably 
generalized - see Q) imparts certain rigid features to the theory of quantum gravity. One basic aspect is the existence 
of quantum observables, which are gauge-invariant, Hermitian operators acting on the Hilbert space of states of the 
theory. An outstanding problem is to understand properties of these observables. 

While a complete discussion of observables in quantum gravity must obviously await more complete understanding 
of its Hilbert space and dynamics, we already have a good deal of information if any such more complete formulation 
is to match onto a quantum version of Einstein’s theory in weak field regimes. One can study the properties of 
quantum observables in these regimes. Moreover, such properties are likely to provide further information about 
the mathematical structure of the theory in the strong-gravity regime. One of the challenges of quantum gravity is 
specifically to formulate gauge-invariant observables that reduce to the usual observables of local quantum field theory 
(LQFT) in the weak-field limit. 

In order to define gauge-invariant observables that reduce to those of LQFT, ordinarily a relational approach is 
taken @-[ 13 , where for example a particle or field operator is localized with respect to some features of the state, or 
with respect to another particle or field operator. Some examples of such constructions are given in We also expect 
that there should be observables that act on a state of the system, say the vacuum, and create or annihilate a particle, 
as in LQFT. These are expected to be the simplest operators reducing to simple operators of a non-gravitational 
LQFT. A key point, however, is that such observables must also create the gravitational field of the particle, in order 
to be gauge invariant (that is, satisfy the constraints). Such operators have been constructed in gauge theories - going 
back to the work of Dirac[l^ - but not, to our knowledge, in gravity. 

Indeed, one way to think of constructing such operators has a close parallel to construction of other relational 
operators 0; we can demand that the point at which a field operator acts is a fixed geodesic distance from a fixed 
feature (or “platform”), which we may take to approach infinity. Such a specification will be diffeomorphism invariant 
for diffeomorphisms vanishing at infinity. An example, in anti de Sitter (AdS) space, is to base such coordinates 
on spatial infinity, which there serves as the “platform.” An equivalent way to describe these constructions defines 
operators by working in a specific gauge, e.g. u sing Gaussian normal coordinates based on the asymptotic platform. 
Such a construction was considered for AdS in |l6l4l^ : a related construction appears in [l9j| . based on earlier work 

[23. 

These observables, acting on the vacuum, create both the quantum associated with the field operator, and a non¬ 
trivial gravitational field. The field for the “Wilson-line” observables we have just described is a singular gravitational 
string. Such a string, once created, is expected to decay to a more natural, less singular gravitational configuration. 
An approach to deriving the operator that directly creates such a configuration is to average over the directions of 
the gravitational string. Working in the linear theory, we will find that such a procedure indeed produces an operator 
that creates the gravitational analog of the Goulomb field, namely the linearized Schwarzschild solution; a parallel 
construction works for quantum electrodynamics (QED), producing the Dirac dressing. These operators, and other 
operators that we construct taking into account possible motion of the particles associated with the field operator, 
give simple example of diffeomorphism-invariant observables in quantum gravity, which we explicitly construct to 
leading order in an expansion in Newton’s constant G. 

Another key question regarding gravitational observables is the algebra that they obey. In LQFT, the field algebra 
closely mirrors the underlying manifold structure[2l| and provides a precise characterization of locality, through 
commutativity of observables associated with spacelike-separated regions. An important question for quantum gravity, 
which appears critical to inferring its underlying structure and its interplay with locality, is that of determining 
the structure of the algebra obeyed by its observables . Since we are able to construct such gauge invariant 
observables to leading order in G, we are able to infer the leading-order structure of this algebra, and we do so by 
explicitly calculating commutators of the operators we have just described. We find that these operators do not have 
local commutation relations,^ due to the gravitational dressings that we have described, and specifically they have 
significant departures from the commutators of LQFT in regions previously characterized by the locality bounds of 
|23l - l25l | . These basic algebraic properties of the theory appear to be the universal weak-field behavior of any quantum 
gravity theory that matches Einstein’s in the weak field limit, and thus should furnish important information about 
the more complete theory of quantum gravity. 

In outline, the next section reviews and extends the discussion of gauge-invariant observables in QED, to set the 
stage for gravity. Section HI describes construction of both “Wilson-line” and “Goulomb” diffeomorphism-invariant 


Here, we differ from claims of |18|l . but for reasons that can and will be explained in the analogous and simpler case of QED. 
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observables in gravity, as well as a generalization valid for moving particles, and performs simple checks that these 
operators create the physically-expected gravitational fields. Section IV turns to the question of the algebraic structure, 
first in QED, where we derive nonzero commutators for various dressings and explicitly reconcile a conflict with [^ . 
by also working out commutators using Dirac brackets. We then exhibit the non-trivial and nonlocal commutators of 
various of the gravitational field dressings. Section V provides conclusions and discussion of further directions. There 
are also three appendices: two with basic formulas for quantization of QED and of linearized gravity, and one with 
other useful formulas. 


II. GAUGE-INVARIANT OBSERVABLES FOR QED 


We begin by reviewing and extending discussions of gauge-invariant observables in QED, which we will take to be 
coupled to a scalar 4> with charge q, 

/:qed = , (1) 

where — iqA^. Gauge transformations act as 

A^,{x) -)> A^,{x) - df,K{x) ( 2 ) 

4>{x) —>■ ( 3 ) 

where the gauge transformation parameter A(a:) vanishes at infinity. The parameter a is a gauge-fixing (more precisely, 
breaking) parameter; special choices are a = 0, Lorenz or Landau gauge, a = 1, “Feynman gauge,” and a = oo, 
restoring gauge symmetry. Further conventions and useful formulas for QED appear in Appendix 

The naive expectation that ^ acts on the vacuum to create particles is confounded by the fact that ^ is not gauge 
invariant. However, one may “dress” 0 to make a gauge invariant, if one defines 

^{x) = V {x)(j){x) (4) 

where the dressing V transforms as V{x) —>• V{x). Following [isl such V’s can be found in the form 


V{x) = exp J d'^x' f^{x,x')Afj_{x') 
Under a gauge transformation 0, this becomes 

r(.) ^ r(.)exp = nx)exp 

and will transform as needed provided 


( 5 ) 

( 6 ) 


d'^nx,x')=5\x-x') . (7) 

Here we use the requirement that gauge transformations vanish at infinity. Note also that in the quantum theory, one 
must give a careful definition of the operator <I>, accounting for possible ordering ambiguities. 

There is a lot of freedom in choosing the function /^. However not all of this is physical freedom. Two functions 
are the same if they agree when integrating against all solutions of the equations of motion. We can think of 
inequivalent dressings, following [^ . as corresponding to different “soft photon clouds” surrounding each particle. 
However they need not be composed only of soft modes; for example string-like dressings, which we are about to 
consider, have divergent energy density. 


A. Faraday, or Wilson, line dressing 

Indeed, a particularly simple dressing is that of a spatial Wilson line, or Faraday line, running to the boundary, 
e.g. along the z direction: 


^W^{x) = (j,{x)e^<l fo°° dsA^i^+sz) 


(8) 
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where z is the unit vector in the z direction. Gauge invariance of an operator may be alternately stated as the 
requirement that the operator creates states satisfying the Gauss’ law constraints, and this operator does so by 
creating an electric field that is localized to an infinitesimally-thin string extending to infinity along the z direction. 

To see that ([5]) creates a line of electric flux, we can use the equal-time commutation relation (see appendix 
[E'‘{x), Ajitj)] = i6jSA')(^x — y) to find [l^ a.t t = t' 

[E^{x),<^wAx')] = -q6‘^ix± - x'j_)0{z - z')<^wAx') (9) 

where xx denotes the components orthogonal to the z direction, xx = (i, a;, y). Note that E^ is negative for positive q 
because cj) creates an antiparticle of charge —q. The other components of E commute with $£). The interpretation of 
this equation is that acting with the operator $ d increases the value of the electric field by the addition of an single 
line of electric flux. 

More explicitly, to find the field created by a general dressed operator $, consider 

A^(a;)$(x')| 0 ) = [A^ix), 5>(a:')]|0) -b $(x')^m( 2 :)| 0 ) ■ ( 10 ) 

If the commutator is proportional to $, this shows that the field after acting by $ differs by A^(a:), given by 

= [A^(x),$(x')] , (11) 


from that of the vacuum. 

To find the future evolution of the field associated to the electric string, we evaluate the commutator to leading 
order in q in the region that is causally separated from the string, but spacelike to the point x'. In Feynman gauge, 
this is given by 


Az{x) =iq dslA^^x), Az{x' -b sz)] = 


qeit-t') 


2TTyJ-{xx - x'j_)'^' 


where e{x) = sign(a;). This corresponds to an electric field which at time t ^ t' is given by 

q\t-t'\ 


Ez{x) = -doA^ = 


27r [-(xx - 


( 12 ) 


(13) 


This solution describes the field lines of the initial string spreading out. There is also a corresponding magnetic field, 
since the electric field is changing with time. 


B. Coulomb, or Dirac, dressing 


The string-like configuration considered above is of course a rather unusual electric field for a charged particle; the 
operator creates a singular electric field concentrated in a single line. We expect such an operator to receive large 
corrections in perturbation theory, and this raises serious doubts as to whether such an operator can be rigorously 
defined in the full quantum theory. One approach to avoiding this would be to thicken out the field into a finite 
diameter tube, with a finite stress energy tensor. A related approach, to which we now turn, is to distribute the field 
lines even more widely - providing a more physical dressing for a charged particle. 

In particular, a Coulomb-like field is anticipated to be more typical, and less singular, configuration for the electro¬ 
magnetic field of a charged particle. Such a field can be found by averaging ([5]) over different directions. Specifically, 
taking x = 0 , we can average 


dzAz —^ 


47r 


d^rtdrEAj = / d'^x 


Ar 

47rx2 


which generalizes for x 7 ^ 0 to give 


^d{x) = (l){t,x)exTp^iq f d^x' x)VD{t, x) . 


47 r|x' — xP 

This is the Dirac [l5l| dressing, which is sometimes rewritten after an integration by parts, 


(14) 


(15) 


d^: 


’ 4 ^ = / ■ 


( 16 ) 
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Using the canonical commutators, one may check that this creates a Coulomb field at time t: 


\ t=t' 


q {x-x'Y 
47r jx — x'\^ 


^d{x) . 


(17) 


For earlier or later times we can use the fact that the commutator \A^{x)^ satisfies the Heisenberg equations 

of motion, with initial data given by the equal-time commutators. Outside the lightcone of x', these are, working to 
linear order in q, the source-free Feynman-gauge equations $£i(a;')] = 0. The gauge-dependent term d^Afj_ 

does not appear in this equation, since d^Afj_ generates gauge transformations (see Appendix!^ and so 

[d^Af,{x), $d(x')] = 0. (18) 


This also means the solution A satisfies the Lorenz gauge condition d^Afj_ = 0. The free held equation with these 
initial conditions has solution 

dTT \X'-X\ 

giving the static electric held of a particle at x'. 

The held conhguration (ITO)) may be continued into the future lightcone of x', but then does not in general satisfy 
the equations of motion. The reason is that the solution there depends on the state of motion of the particle at and 
after t'. In general, this solution will have, in addition to (IT^ . a radiative part. The full determination of this held 
depends on particle motion in the interacting theory. 

As we have noted, if the Wilson line operator acts on the vacuum, the highly-localized electric held it creates 
is not stable, and will decay Bin. Another way to see its relation to the Dirac held is to rewrite the held arising 
from the Wilson line as a longitudinal, Coulomb, piece, plus a transverse, radiative, piece [s^. using 

A, = Ai + Aj = ?^A, + (i„ (20) 


The integral f dsA^ then gives the Dirac dressing (IT6l) : f dsA^ contributes extra held energy which radiates to 
inhnity. 


C. Worldline dressing 


An electromagnetic dressing that takes particle motion into account is dehned by the world-line expression 


^*^^( 3 ;) = UD(a, 0)f'exp|-ig J dA^A^(a:''(A))| ^(x) . 


( 21 ) 


Here Vd is the Dirac dressing given in ([15]), T anti-orders in time (latest time is to the right), and x^{X) is a trajectory 
with x^{0) = (a,0) and a:'^(l) = x^. This produces the electromagnetic held that results from classical motion along 
a:(A). For large m, a quantum particle approximately follows such a trajectory, and (I2ip produces the leading-order 
held, but for general m there are quantum corrections to this held. 

To see that (pil) produces the held of the moving particle, we consider again the commutator A^. This will contain 
a term from Vd as well as one from the worldline, given by 


f dx’" 

A^{x'),-iq dX^A,ixP{X)) 


dXD^,{x' -x{X))'^+0{q^) , 


( 22 ) 


where is the Pauli-Jordan function for A^, (IA22I) . When x' is spacelike to x and to the future of (a, 0), D^i, 
agrees with the retarded propagator so this term produces the held of the moving particle. When x' is spacelike to x 
and to (a, (1), we have the Coulomb held (HSj). 

We note that although the vector potential A^ is discontinuous across the lightcone of (a,()), this is a pure gauge 
discontinuity and will be continuous provided the trajectory a;(A) is initially at rest. If A^ is initially at rest, 
then the held just to the future of the lightcone of (a, 0) is 

g 

dTrjx — x'\ ’ 




( 23 ) 
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which is related to the configuration (1191) (with x x') by a gauge transformation 


A(x') 


q{t -1') 
47 r|a; — x'\ 


(24) 


Let us examine the large-m case more closely. In the large-m limit, the four-velocity of a particle is constant, 
and superselectedfs^ . Then, from the equation of motion 




(25) 


one may show 


iu^^d^<^WL{x) = m^wL - ■ (26) 

Thus, to leading order in l/m, ^wl creates a state that is an energy eigenstate; this statement was called the “dressing 
equation” in 

The latter feature, together with the fact that in the asymptotic regime, where particles have large separation, m 
is the largest energy scale, motivated [13,113 to use m, with constant-velocity trajectories, to define asymptotic 
states which avoid the usual infrared difficulties of describing asymptotic particle states in QED. 


D. Lorenz dressing 

One natural way to specify a gauge-invariant operator is to fix a gauge by some condition, and define $(a;) to be 
(j){x) in that particular gauge. If this condition completely fixes the gauge, we can solve for the gauge transformation A 
that transforms an arbitrary vector potential into this chosen gauge. This allows us to express $(a;) in a manifestly 
gauge-invariant way, 

$(x) = (27) 

where X{x) is a nonlocal function of Indeed, the preceding construction of for example, gives the field 

operator in axial gauge. 

Another apparently natural choice for this purpose is Lorenz gauge = 0. However this is not a complete 

specification of the gauge: given a vector potential we must find the gauge transformation A to Lorenz gauge, 
which is found by solving 

d^^A^ = DA. (28) 

The solutions to this equation are not unique until we specify appropriate initial (or final) data for A. 

If our initial data surface is spacelike, we should choose initial data that determine both A and its first derivative 
away from the surface. For example, we can fix n^A^ = 0, where is a normal to an initial data surface, which 
we take to be a spacelike of constant time. This leads to an initial condition for n^9^A = This does not yet 

determine a solution for A, since we need to know both the initial value of A and its time derivative. We can give 
an additional constraint = 0, which leads to the equation V^A = d^Ai on the initial surface. This determines A 
completely, once we include the boundary condition A —>■ 0 at spacelike infinity. 

Of course there are other ways to fix the gauge asymptotically. Since A satisfies a free wave equation, it may make 
more sense to fix its initial data on J^“. For example, we can impose Au = 0, where u is the null generator of 
This leads to the equation 9„A = A^, which then determines A completely on from the condition that the gauge 
transformations vanish at infinity, A = 0 on 

One can therefore likewise formulate dressings based on these gauges; we leave description of these dressings to 
future work. 


III. GAUGE-INVARIANT OBSERVABLES FOR GRAVITY 

We next turn to a discussion of gravitational observables analogous to those we have described in QED. The 
nonlinearity of gravity of course is a significant complication, which we manage by working perturbatively in Newton’s 
constant Gu. In particular, the leading-order, linear structure of gravity has important similarities to that of QED. 
This structure will describe the dominant behavior in the weak-field regime. 



( 29 ) 


Specifically, consider the gravitational Lagrangian, minimally coupled to a scalar (j) with mass m, 


r 2 1 VM 1 

^ 1^1 


v° 




where = 32ttGd, and R is the scalar curvature. The derivative V*’ is with respect to the background metric about 
which we perturb. The parameter a is a gauge fixing (more precisely, breaking) parameter; a = 0 gives de Bonder 
gauge, an analog of Lorenz or Landau gauge, a = 1 is an analog of Feynman gauge for QED, and a —>■ oo is the 
analog of unitary gauge, which restores the gauge symmetry. The perturbative expansion follows from 


= 




(30) 


where here the background is taken to be the Minkowski metric. Further details of the perturbative theory are 
supplied in appendix [B] 

As was the case for a charged field in QED, the field (j){x) is not a gauge-invariant operator, since under a diffeo- 
morphism f M —>■ M, the field transforms as the pushforward 

{f*4>){x) = 4>{f~^{x)) (31) 

If we take / to be an infinitesimal diffeomorphism, f^{x) = 

5(j){x) = (j){f~^{x)) - (j){x) = + 0 {k^). (32) 


To find a diffeomorphism-invariant operator, we will form a composite operator that includes both <f) and the gravi¬ 
tational field sourced by (j). 

A simple approach to constructing operators invariant under linearized diffeomorphisms, analogous to the approach 
taken for QED, is to seek an operator of the form 

<I’(a;) = = (j){x^ + V^{x)) = (j){x) + V^{x)dfj,(j){x) + 0{V'^) , (33) 

where (compare the QED expression ([3|)). Here the “dressing” V^{x) is a functional of the metric 

perturbation that transforms under a diffeomorphism as 


x^ + V^{x) —>■ f^{x + V{x)) , 


which at linear order becomes 


(34) 


SV^^ix) = <^(a:) . (35) 

Then, using the combined transformations, at linear order 

5$(x) = S[(l){x + H)] ~ 6(l){x) + SV^{x)dfj,(j){x) = 0 . (36) 

The transformation law (1351) should follow from the change in the metric, which transforms as 

~ ~ “b ^J 

correspondingly the metric perturbation transforms as 

+ 0{k) . (38) 

This paper will primarily (though not exclusively) consider such constructions at leading order in k, and in so doing 
will consider only linearized diffeomorphisms and drop the 0 (k) term from (1381) . A construction to higher-order in k 
will then be constrained by this leading behavior; we save general examination of all-orders behavior for future work. 
As in the case of gauge theory, we can try to find E’s satisfying (1551) that are of the form 

V^{x) — ^ J f^''^{x,x')h^x{x') , (39) 

where / is assumed symmetric in its last two indices . Given (1551) . in order for to transform as dSSl), 

the function / parameterizing the dressing must satisfy 

2d'Jf^''^{x, x') = 5'^{x - . (40) 

Again, there is substantial freedom in choosing the functions /, which determine the “soft graviton cloud” (plus, 
possibly, a harder component) oi a (j) particle, though not all freedom yields physically-inequivalent dressings. 
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A. Gravitational Wilson line 

As in QED, a very simple dressing is a Wilson line. One very geometrical way to think of such a line in gravity arises 
if points are located by shooting a spatial geodesic in a perpendicular direction from a fixed asymptotic “platform,” 
where diffeomorphisms are taken to vanish; ultimately, we might take this to reside at infinity.^ Specifically, the 
location of a general point is determined by the platform point at which the geodesic originates, and the distance 
along the geodesic. The field (j>, expressed as a function of this data, is gauge invariant under gauge transformations 
vanishing at the platform. 

We can thus choose as coordinates the initial platform position, and the geodesic distance. Specifically, define new 
coordinates in which the z direction is perpendicular to the platform, and choose the “axial” gauge (or, geodesic 
normal coordinates) 


h,f,=0. (41) 

In this gauge, the z coordinate is the geodesic distance, and the remaining coordinates correspond to the platform 
position. The diffeomorphism-invariant held is thus just the scalar held (j) in these coordinates. To write this in terms 
of the scalar held in arbitrary coordinates, consider a new (general) coordinate system 

x^ = r{x) . (42) 

Note that, for general coordinates x^, the function depends nonlocally on the metric, since it must take a general 
metric to the form (1411) . Then, 


(43) 

is the diffeomorphism-invariant held written in terms of the held (j) expressed in general coordinates. (At this stage 
we drop the accent on x.) 

The preceding statements are valid to all orders in the perturbative expansion in k, but can also be simply explained 
to leading order in this expansion. Specihcally, consider a small metric perturbation and parameterize the relation 
between axial coordinates x^ and a general coordinate x^ as 


Then, the relation (1431) becomes 


x'^ = x' 


-VA 


Wz 


(44) 


4>w, (a;) (x)) (45) 

to leading order. Given the linearized gauge transformation (I38L the axial gauge (ED can be hxed by a diffeomorphism 
of the form 


= 2 y dsh:,^{x+sz) = Vw,,z{x) 



hzpiA + sA +2^1^ 



ds' hzz{x + s'z) 


(46) 


where jl denotes indices excluding z, and the platform has been taken to z = oo. 

With these expressions, one can explicitly check the leading-order diffeomorphism invariance of (ED, (ED- From 
the metric transformation (I38p , transforms as in (1351) for a diffeomorphism ^ that vanishes at inhnity, so according 
to (1551) . is diffeomorphism invariant. The gauge-invariant operator is the gravitational analog of the electric 
Wilson line operator dehned in ([HD. 

This construction can also be characterized by solving the geodesic equation for a curve x*^ = x^ -|- sz^ + v^{s) with 
the boundary condition v^{oo) = 0. The quantity V^{x) = 'u^(O) then gives the coordinate displacement between 
general coordinates and axial coordinates. To find the linearized dressing , we solve the linearized geodesic equation: 


dlv'^is) + szA = 0 


(47) 


which gives 


= - dssrA{x + sz) + s.t. , 


( 48 ) 


^ For an analogous construction in anti de Sitter space, with the AdS boundary providing a platform, see [T3 . Il8l . [3^ . 
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where s.t. denotes a surface term at infinity.^ This can be checked to be equivalent to (1461) . 

As with the case of electromagnetism, the operator (H51) . (1451) creates a singular metric configuration, where grav¬ 
itational field lines are concentrated in an infinitesimal string. This can be seen by studying the commutator of h^i, 
with . Even in the linear theory, this is an unstable field configuration, and if it is created the field lines will then 
dynamically spread out into a Coulomb-like configuration. Of course in the non-linear theory we also expect such a 
concentrated gravitational field to source large non-linear corrections, and as a result at a minimum the gravitational 
field configuration should be “thickenedjl^;” the absence of string-like solutions to the sourceless Einstein equations 
moreover likewise strongly indicates that the field configuration should undergo such a decay to a more symmetric 
configuration in the full theory. 


B. Gravitational Coulomb dressing 

1. Construction 


Like with QED, an approach to finding a more symmetric dressing is to average the gravitational Wilson line over 
all directions. This is most easily done starting with (I48L which then becomes, specializing to a; = 0, and working in 
D = 4, 


^c(0) = -^ J dr 




I Kpf-fP+s.t. 


(49) 


with f the unit radial vector. While we use D = 4 for much of the following, the discussion readily generalizes to 
D > 4. Specializing to the timelike component, we find 




where the surface term has cancelled. The spatial components become 

where again the surface terms cancel. At a general location x, these become 




K f ,3 , [ d‘^d^dohap hoad° 


2d 


d? 


, Vh[x) = 


3k 

Stt 


^3 




^OLjS 


(50) 


(51) 


(52) 


where d = x' — x. One can check directly that Vq has the correct gauge variation; using (1381) , one finds 

3k C ^ A A 

SV^(x) = kC{x) d^^'-^f' -dd-id^ 


f d^n'r'^f' -d 


5Vq{x) = k5°(x) -I- 


47r 




(53) 

(54) 


Thus we recover the expected transformation law, 5Vq{x) = k^^{x), as in (l35l) . for ^ satisfying falloff conditions such 
that the boundary terms in these expressions vanish, i.e. and r^r + ■Co vanish as r —>■ oo. 


2. Dressing field 

Comparing QED, we might expect that the operator 4>c(a;) = (/>(x -I- Vc{x)) creates a gravitational analog of the 
Coulomb field. In order to check this, we examine the commutator [h^jy(a:), 4>c(0)], which indicates, as with the QED 
case, how the gravitational field is changed by action of the dressed operator. 


® This vanishes if \x\hij,z{x) vanishes in the limit x oo. 
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As a first step, we need the commutator 




47r 




t f 


IK 

47r 


2r 


= / d\ ( 




] 11(2;' - x) , 


(55) 


where we work in Feynman gauge and D is the Pauli-Jordan propagator of a massless scalar, (IA20I) . This can be 
rewritten in terms of scalar integrals using (see appendix ICj) 




- dad^r, 


pi 

9 — 

r>^ 


(56) 


where the spatial metric is denoted We then integrate by parts, so that the derivatives act on D{x' — x), and 
then trade them for x' derivatives, which can be pulled outside the integral; 


IK 


[hf,^{x'), V'^(O)] = ^ I d^x[ ^d^d^rdo - ^do + ) D{x' - x) 


47r ^ 

p/ p/ 0/ 
47r M t 


d^3 




D{x' — x) 


(57) 


where d denotes the spatial gradient. The resulting scalar integrals can then be evaluated: 


/ 


-1 ifF>r' 

d^x - D(x' — a;) = ^ —t'Jr' if — r' <t' < r' 
1 if t' < -r' 


(58) 


jS.rDW-.) 


r-t'2-2^ iff'>r' 

+ ift'<-r'. 


(59) 


Using these result with (1571) and combining terms we find, in the region that is spacelike to the origin, —r < t < r, 


[h^,{x),vsm 


Ait 2r 2r^ 


(60) 


In the regions that are timelike to the origin, t > r and t < —r, we find [hfj_^{x), U(9(0)] = 0. We can check explicitly 
that (1501) satishes the harmonic gauge condition = 0 as well as the equation of motion = 0. 

We also need the commutator [/i^i^(a;), U(3(0)]. Only the spatial components contribute, and they are given by 

[hki.x'),V^{Q)\ = ^ J - x). (61) 

We use the same trick of expressing the tensor as a derivative (see |c3), 

= ^d^djdkr - q^.jdk) i , (62) 

integrate by parts, pull the derivatives outside the integral, and find 

— *^7 K / "t 1 \ 

[hMx),Vim = (63) 

in the case where x is spacelike to the origin. In the timelike regions t > r and t < —r we have [hfj,^{x), U(1(0)] = 0. 
We can again check that the result (100)) satishes the harmonic gauge condition = 0 and the equation of motion 

= 0 . 
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We first consider the dressing field due to a massive particle, which we can approximately consider to be at rest, so 
that do(j}{x) = im(f>{x) (for the creation part of the operator) and we can neglect spatial gradients. Then we have 


[hf,^{x),^{0)] = [hf,^{x),VS{0)]do(l){0) = 


47r 


2r 2r^ 




^(0) = h^j.v{x)(j){Q) 


(64) 


which gives the field of the particle, valid for all points spacelike to the origin, |t| < r. As noted, this satisfies the 
harmonic gauge condition and the equation of motion \I\hfii, = 0 in this region. 

While the metric in (1641) may be unfamiliar, it is simply the linearized Schwarzschild metric in an unusual gauge. 
By means of an infinitesimal diffeomorphism 


rx L , ■] o / 




(65) 


we can put the metric into the form: 




Km 




which is the linearized Schwarzschild solution. This can also be checked by calculating the ADM energy, 

1 


pO _ 
^ADM — 


IbTrG 


j) f^{djgij — diQjj) ~ ~ r^{djhij — dihjj) . 


( 66 ) 


(67) 


which indeed yields Tadm = the solution in (l64ll . 

Generalizing to the case of a localized source with d^cj) ~ we also pick up a field contribution from ()63p . of 

the form 


3k 


Stt 




1 


= -ttP" + -d^djdk- 


9 


( 68 ) 


Again, this satisfies the harmonic gauge condition and the equation of motion = 0 in the region spacelike to the 

origin. 

As a check, we can calculate the ADM momentum of the solution. This is given by [3^ 

C^DM — ^ T d ^ G^ADM — ^ T d ^ ^j^Oi ; (h9) 


where the linearized canonical momentum in the ADM formalism: 

^ADM — ^i^Oj djhoi hj^]^6ij T 2jdkhQi^Sij , (^t)) 

and we have used the identity f r'^d^fl {Okhok^ij) = f r'^d^fl {dihoj). We can verify that the solution (155)) has 
doM =7**1 S'® expected. This is a general consequence of the commutation relations 

• (71) 

Another check on the preceding constructions is to consider the commutator 

Kdm: = id>^<^{x) + 0(k) . (72) 

This shows that the ADM momentum generates translations of the diffeomorphism-invariant observables. Note, 
however, that at this order Padm only generates translations of the field tp and not of the dressing; this is a consequence 
of truncating our perturbation theory at 0(k). Since Padm contains an explicit factor of k~^, one has to go to order 
to see the effect of translating V in the operator $ at order k". We expect that the ADM D-momentum will 
continue to generate translations of $ at higher orders in perturbation theory. 
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3. Further comments on different dressings 


We have noted above that the gravitational Wilson line will in general decay to a Coulomb-like configuration. 
As with QED, the averaging procedure yielding projects the Wilson-line field onto this configuration, and the 
remaining transverse field, as in (1201) . represents the radiative part of the gravitational field. 

Other simple examples of different dressings satisfying the constraints can be obtained by adding to a giving a 
solution of the constraints a multiple of the deformations 


Stt , 


AC°(a;) = -3— / d'^x' -AJl— 


(73) 


or 






T-jk 


(74) 


again with d = — x. These can be shown to be invariant under diffeomorphisms (|38p which vanish sufficiently 

rapidly at infinity, and thus adding them maintains the condition (1351) : moreover they shift the dressing field by a 
diffeomorphism. Note that a linear combination of them can be added to the Coulomb dressin g to cancel the leading 
f dependence at infinity, resulting in invariance under supertranslations at spatial infinity [381 l39j . 


C. Gravitational worldline dressing 


In parallel with QED, the dressing creates the field for a particle at x that has been at rest forever. Note 
also that it does not necessarily give the correct field in the future light cone of the point x, since that will depend 
on the subsequent state of motion of the particle. We can consider gravitational dressings corresponding to more 
general trajectories for the particle; for example, for an operator describing a particle at x, we might want to consider 
gravitational dressings corresponding to different paths that the particle took to x. For completeness, we give a 
construction relevant to these cases. A possible way to account for this is to construct a dressing similar to that of 
m, beginning with the Coulomb dressing (j49|) at a time in the distant past, and then adding a worldline component 
to account for the motion of the particle up to a given time. In principle we might try to consider an arbitrary 
worldline, specified in an invariant way, e.g. by specifying its acceleration in a local frame carried by parallel transport. 
However, this will not yield a dressing creating a field solely generated by a particle following that worldline, since 
the gravitational field must be coupled to a conserved stress tensor, which for an isolated particle must correspond 
to a geodesic. Thus allowed worldlines are determined by the geometry, and the final position and momentum of the 
particle. 

Alternatively, we can proceed directly to construction of the worldline-dressed operator directly analogous to (1451) , 
via a geodesic construction. Specifically, we define a worldline-dressed operator of the form 

^wl{x) « cl){x^ + V^^{x,x') + V^{x')) (75) 

where the ~ denotes that we are again working only to linear order in the metric perturbation. Here Vq{x') dresses 
the point x' = (t',x) directly to the past of x = (t,x), and transforms as 5Vq{x') = k^^{x'), as before. This is the 
special case of a geodesic with zero initial velocity, and will be generalized below. Then the worldline part of the 
dressing can be associated with the choice of a gaussian-normal gauge = 0 for t > to, analogous to the axial gauge 
choice. Effectively, one localizes spatial points at t = t' with respect to the boundary, and then localizes future points 
by relative to these. As in eq. (l46l) . this gives an expression 


K f*~*' 

VwlA^,x') = / dXhooix' + Xi) 

^ Jo 


0 

t-t' 


VwL,iix,x') = - f dX KhQi{x'+ xi) - l-di f dX'Khoo{x'+ X'i) + diVc,o{x') 
Jo ^ Jo 


(76) 


The latter term in is needed so that 


dV^j^{x,x') = <'"(x) - <^(x') 


( 77 ) 
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as is readily checked. Then, under a general diffeoinorphism 

5<^wl{x) « + V^^[x,x') + V^{x')) + d^c^{x^ + V{^^{x,x') + V^{x')){5V^^ + <5^^) = 0 , (78) 

as in (|36ll . 

These expressions can be derived directly from a geometrical construction, while simultaneously generalizing to 
non-zero initial velocity. Specifically, we first use the Coulomb dressing to establish a frame at the point x' = {t',x') 
on the surface defined to linear order by 

x>^=x>^ + V^{x^) (79) 

with constant t' and varying x*'. The spatial vectors of this frame are given by (dropping accents) 

e'^ = S^ + d.VS. (80) 

The timelike vector Cq is giving by finding the unit normal to the surface defined by (179^ : the conditions 

0 = +'^hfj,^)e^di{x''+ Vc) = KhQi +eoi +diVco , -I = (81) 


give 

^0 = (1 + + 5^\-K,hoi - diVco) ■ (82) 

A geodesic is determined by shooting it from x' = (t' ,x') with an initial four-velocity specified with respect to this 
frame. The choice in (1751) corresponds to an initial four-velocity purely in the normal direction, doV^j^ = Cq. More 
generally, we can take initial velocity 


drV^{x, x') 


(83) 


Solving the geodesic equation 

+ ro, = 0 (84) 

for the perturbation from the straight-line trajectory, with these initial conditions, gives the curve x^ = x^ + 
y^^(x,x') -I- Vp(x'). This determines the worldline dressing 

y^^(x,x')=/ dA[(A-r)rO„(x'-f AM)-fu''eO(a;')] . (85) 

This can be checked to be gauge invariant, beginning with the gauge transformation for the frame field: 


(5ef = d,{5V^) = 

= -K(do^o)t^ + - udi^o) = . ( 86 ) 

Then, (1851) can be varied; a u-dependent term from varying the Christoffel symbol is cancelled by the second term, 
using (1551) . and we again find the necessary variation, (1771) . 

In summary the combined Coulomb and worldline dressings can be understood by the following procedure. First 
we establish an equal-time surface at time t' by the geodesic averaging/Coulomb construction, then to locate a point 
at a later time t we shoot a geodesic forward from this surface for a given proper time in the direction specified 
with respect to the frame at the surface. One may alternately add a worldline dressing to the Wilson line-dressed 
operator, analogous to (ITSl) . 

The leading-order dressing field resulting from the worldline construction can also be worked out, in analogy to the 
calculations of IIII B 2l We leave the description of this for future work. 


D. Gravitational dressing equation 

In Ref. [3^, dressed field operators for QED were derived from the requirement that in the infinite-mass limit, the 
equation of motion for the dressed field should reduce to a first-order equation This first-order dressing 

equation can be shown from the equations of motion to be solved by the worldline dressing m, up to terms of order 
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1/m - see (l26l) . The construction depends on a choice of Lorentz frame, so it selects a preferred dressing for each 
frame, corresponding to different boosted Coulomb fields. 

Also for completeness, we here show that a similar equation is satisfied by the gravitational worldline dressing. If 
we study the large-m limit, that means we consider solutions with a rapidly oscillating phase where t will be 

the proper time coordinate used in defining the worldline dressing for gravity. Here we work in the rest frame of the 
massive particle, corresponding to choosing a geodesic with zero initial velocity in the preceding subsection, although 
the discussion may be generalized to non-zero velocity. The equation of motion for (f> can be rewritten as 

»■ p-imt 

- = m^- — -V2(e™V)- (87) 

2 2m 

Here we have used the fact that t measures proper time, = — 1. To rewrite this as a gravitational dressing 

equation, we need to reexpress it in terms of the dressed field, and to relate the first term to a time derivative djdt. 
Recall that the worldline dressed operator ^wl{x) is defined as the value of the field at a set of coordinates x with 
an invariant physical meaning. Like with the gravitational Wilson line, t is the proper time along a timelike geodesic, 
and X labels the different geodesics. Thus 


^wl{x{x)) = (j){x) ; (88) 

compare (1421) and (l43ll . Using this, we can reexpress the first term in (l87l) 

- = -i{V^t){Vp,x'^)^^wL{x{x)). (89) 

Since the worldline coordinates are Gaussian-normal, (V^t)(V^x‘") = rf", and the equation of motion (l87l) becomes 

i^^WL - l:{^^i)^WL = m^wL - ^7T^V2(e™‘$rvL) , (90) 

at 2 2m 

where d/dt is defined with the coordinates x* held fixed. This is first-order in time, up to 0{\/m) corrections. 
Dropping the C>(l/m) term gives the gravitational dressing equation, analogous to that discussed for QED in [s^. 


IV. ALGEBRAIC STRUCTURE 

An important basic question for a theory is the algebraic structure of its algebra of observables. In local theories 
without long-range fields, the algebra has a net structure of subalgebras that mirrors the decomposition of the 
underlying manifold into spacetime regions. Long range gauge fields lead to additional subtleties with such subalgebras. 
We next turn to an examination of some basic aspects of such algebraic structure, focussing on the gauge-invariant 
observables that we have constructed above, beginning first with QED and then turning to gravity. 


A. Algebra for observables in QED 

1. Commutators for Faraday dressing 


Let us first consider the scalar field operator dressed by a Earaday line introduced in section III Al If we consider 
commutators of (x) with (x') at equal times t = t', they will vanish because the operators only involve and 
Az, all of which commute at equal times (see appendix |A] for basic commutators). However, away from equal times 
there will be a nonzero commutator whenever any part of the two electric strings are causally separated. 

To show this we consider one point slightly to the future of the other, which is determined by {x),^w^ (a;')] 

at equal times. Using the definition of (|8]), and the identity 


dt 


e°=e° 


o + i [d,o 


whenever [O, O] is a c-number, we find that 


^wAx) = 


^{x) + iq(j){x) / dsAz{x -\- sz) + -zA4’{x)5^A)) / ds 
Jo 2 Jo . 


(91) 


(92) 
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The last term arises from the commutator (see eq. (IA21I) 1 Az{x')] = —iS^{x — x'); it is formally infinite, and can 

be interpreted as due to the infinite energy of the singular string. It can be regulated by cutting off the construction 
at z = Z, and by smearing out the singular string over a small region. However, since it is proportional to ^ it can 
be alternately absorbed by a c-number phase rotation of ^ and thus also does not contribute to the commutators 
of present interest. 

If we now commute (j92p with <I>w^(x') for x ^ x', the contribution of the first term vanishes by [(/>, (j)] = 0, and we 
are left with 


[^wAx),^wAx')\ = ds[AAx +sz),^wAx')] 

Jo 

noo 

= iq'^^Wz{x)^w^{x')S^^Ax± — x'j_) / dz”. (93) 

J m.six.{z,z') 

The last term is infrared divergent and may be again regulated by cutting off the z integration at finite Z. 

To understand the appearance of the nontrivial commutator, imagine first acting on the vacuum with the operator 
^w^{x). This will create a charged particle, and a nontrivial configuration of the electromagnetic field, ©, (TT^ . If we 
subsequently act with at a slightly later point {t + 6t, x), there is a large phase associated with the introduction 
of the new charged particle into the preexisting nontrivial field. The divergence arises since this field does not decay 
as 2 > oo. While the Faraday dressing provides a nice intuitive picture of the origin of the non-zero commutator, the 
divergence can be eliminated by working with dressings like that of Dirac, which are better behaved. 


2. Commutators for Dirac dressing 


Equal-time commutators of Dirac-dressed operators (1151) likewise vanish at equal times, but can be shown not to 
vanish as one operator is moved into the future by considering the commutator [$d(x), ^^(x')]. Specifically, the 
needed time derivative is 


^d{x) = Vd{x) ( Ax) + iqAx) / d^x' 


jx' - xf 
47 r|x' — xP 


AAt,x') + Ttl^Ax) / d^x' 


(47r)2 


(94) 


Here, again, the last piece is interpreted in terms of field energy; here it is infrared finite, but ultraviolet divergent, 
though does not contribute to the commutator of interest. The equal-time commutator arises, as with the Faraday 
dressing, from the second term in (IMl) . and takes the form 


[4>D(x),$D(a;')] 


iq'^^D{x)^D{x') 


ix" - 
47 r|x" 

f dA" 



{x" - 
47 r|x" 


[Hi(t,x"), 4>z>(a:')] 

- xy [x" - x% 
— xP 47r|x" — x'P 


iq^^D{x)^D{x') J d^x”- 


47 r|x — , 


iq^ 


47 r|x — 


— $I5(x)$d(x') . 



(95) 


In contrast to the divergent commutator of the Faraday-dressed operators, the commutator of the Dirac-dressed 
operators is nonsingular and decays with distance, reflecting the fact that the Coulomb field spreads out with distance. 
Moreover, we see that the coefficient entering the commutator is precisely the Coulomb energy between charges at 
points X and x'. 


3. Relation to calculations with Dirac brackets 

Note that claims that commutators of the field </), in axial gauge, vanish at spacelike separation, and argue that 
this means microcausality is preserved. This seems in conflict with the results derived above. Here we will reanalyze 
the question, in the Dirac bracket formalism used in [^ . and show that the commutators in question are indeed 
nonvanishing. To do this we review the Dirac bracket formalism!^, and compare it with the construction used in 
the present paper. 











17 


First, we recall that the dressed operator coincides with the value of (j){x) in axial gauge = 0. So 

the commutation relations of {x) are given by the commutation relations of 4’{x) in axial gauge, which can be 
evaluated using the Dirac brackets. We will encounter a possible ambiguity in this gauge choice, since can really 
only be set to zero by a gauge transformation vanishing at infinity if / dzA^ = 0; correspondingly, the operator 4>{x) 
could be accompanied by an electric string running oS to z = oo, or to 2 = —oo, or a linear combination of the two, 
and these are gauge-inequivalent configurations. 

To carry out the Dirac quantization, we begin by treating the components of A^ as canonical variables, and find 
their momenta tt^; these are given in (IA9I) with a —> oo. Only tt* are nonzero, and = 0 gives a primary constraint. 
In order that this constraint be preserved in time, we find the Gauss’ law constraint, = —j^, where is the 
charge density defined in (I ATI) . These constraints Poisson-commute with each other and with the Hamiltonian, and 
give a system of first-class constraints. 

In order to fix a gauge and reduce the system to physical degrees of freedom, we add an additional gauge-fixing 
constraint, A^ = 0. This constraint does not commute with Gauss’ law, or with the Hamiltonian. To preserve the 
constraint in time, we introduce its time derivative tt^ + dzAo = 0 as a further constraint. Note that this constraint 
does not commute with the primary constraint = 0. We thus have a system of second-class constraints, given by 

7r° = 0 , 9i7r* -I- = 0 , = 0 , + dzA[) = 0 . (96) 


Letting Xi label the constraints (there are 4 per point) we can consider the matrix Cij = {Xi;Xj}; written in terms 
of the Poisson brackets {A^{x),'k'^{ x')} = 5''^5^{x — x'). If this matrix is invertible, then we can form its inverse , 
and define the Dirac brackets. 


{A, B}db = {A, B} - {A, Xi}C^^{x,.B} . 

Because the term C*-! is obtained by inverting a differential operator the Dirac bracket is nonlocal. 

In the case of axial gauge, the nonzero Poisson brackets are 

{7r°(a:),7r^(a;')-I-92Ao(a;')} = dz5^{x-x'), 

{diTT^x) + f{x),Az{x')} = -dzS^{x-x'), 

{Az{x),'k’'{x') - dzAo{x')] = 6^{x-x'). 

The constraint matrix C can be expressed as a 4 x 4 matrix of differential operators, which takes the form 

■ 0 0 Odz- 

0 0 -a, 0 

0 0 / ■ 

.dz 0 -/ 0 . 


(97) 


(98) 

(99) 
( 100 ) 


( 101 ) 


The general antisymmetric inverse of the constraint matrix is given by 


C ^{x,x') = 


0 


i|z - z'\ +cz + dz' + e 

0 


— ^\z — z'\ — cz' — dz — e 


ie(z-z')- 


ie(z - z') +1 


-\e{z-z')-d 


0 

-\e{z-z') + d 0 

0 0 

0 0 


d^(x_L - x(l) 


( 102 ) 

We have to choose boundary conditions in order to invert C, which determines the constants c, d and e. The freedom 
in inverting the matrix of Poisson brackets is due to the ambiguity in transforming to a gauge such that A^ = 0. 
The Dirac bracket defined in corresponds to the choice c = d = e = 0. These different solutions define different 
dressed scalar operators. The simplest way to see this is by considering the Dirac bracket of (j) with the z-component 
of the electric field, which measures the electric flux. 

To find the Dirac bracket of (j) with E^, we need the nonzero Poisson brackets of these quantities with the constraints, 
which are given by: 


{diTT^x) + j°(x), (/)(x')} = iqS^{x - x')cj){x') , 
{E^ix), Azix')} = 5^{x - x') . 


(103) 

(104) 


The Dirac bracket is then given by 


{E^x),(j}{x')}DB = - J d^x"d^x"'{E^x),A^{x")} [-^e{z" - z"') - d] 5‘^{x'[_ - x'Didiir^x'") + f{x"'),(j){x')} 


= iq [^e{z — z') -f d] S'^{x± — x'j_)(j){x') . 


(105) 
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Thus we see for d = 1/2 we have a Faraday line of flux —q pointing in the z direction, which is the dressing of 
© with the identification [ , ] = i{ , }ob- The choice c = 0, used in gives a particle dressed by two Faraday 
lines of flux q/2 pointing in opposite directions. 

We can now consider Dirac brackets between the scalar field (j) and its canonical momentum The modification 
entering the Dirac brackets comes because of the commutators of the matter field with Gauss’ law is nonzero. However 
C'-f is zero when i and j both label the Gauss’ law constraint, so these Dirac brackets coincide with the canonical 
Poisson brackets. This conclusion was reached in (^ . and used as part of an argument for commutativity of (^’s at 
spacelike separation. 

However, while this argument shows that the Poisson brackets between (j) and tt^ are unmodified by the coupling 
to the gauge field; it does not show that </> commutes with ^ at a later time. To see this we can consider the Dirac 
bracket of with (f>. This leads to a nontrivial Dirac bracket due the Poisson bracket of Aq with the 

primary constraint 

{^{x),(I}{x')}db = - J d^x"d^x'"{iqAo{x)(j){x),TT°{x")} - z"'\ - cz'" - dz" - e) - x'”){f{x"'),(/{x')} 

= —q^ (||z — z'\ + cz' + dz + e) 4>{.x)4>{.x') ■ (106) 

This indeed agrees with the commutator of the Wilson line dressing (IMl) . if we take c = d = \ and e = —Z. 

Thus, in summary, while one does find that the equal-time Poisson bracket {4){x),'K(j){x')'\ = 0 lor x ^ cc', this does 
not imply that {4>{x), (j){x')} vanishes at spacelike separation, so the operators do not commute in general outside the 
lightcone. 


B. Algebra for observables in gravity 

We finally turn to the important question of the algebra of gauge-invariant operators in gravity. In section IHII we 
described the construction, at linear order in k, of diffeomorphism-invariant observables corresponding to a matter 
field and its gravitational dressing. In parallel with the preceding discussion on QED, we can now likewise investigate 
the algebra obeyed by these observables. While these have only been constructed to leading order in k, and thus we 
will only find the leading-order commutators, these have interesting structure, and of course constrain the all-orders 
commutators since the latter need to match our results when expanded to this leading order. 

The general form of the gauge-invariant observables we have described is 

$(x) = (j){x) + V''{x)d^4>{x) + K^$( 2 )(a:) -I- , (107) 

where is of order k, and explicit examples have been given in section Hill We will consider the equal-time commu¬ 
tators [<!>,$] and [<!>,$] to leading order, which is to order At this order and for x ^ x' these commutators take 
the form 


[$(a;), $(a;')] = [V^x), V'{x')]df,(j){x)d^(j){x') (108) 

and 

[4>(a:), 5>(a;')] = [V''{x),V''{x')]dfj,(l){x)d^cj){x') + [Vix), V{x')]dfj,<j){x)d^(j){x') . (109) 

Note that validity of (|I08|) and (11091) requires elimination of other terms that could potentially contribute at order 
First, there is a possible contribution to the commutator (llOSp of the form 

[V''{x),cj}{x')]df,(l){x) (110) 

and similarly with x ^ x'. However, V^{x) as we have constructed it depends only on the metric and its first 
derivative on the constant time slice, so this commutator vanishes. 

There is also the possibility of a term in (11091) of the form 

[V''{x),4>{x')]d^4>{x), (111) 

and likewise with x x'. Indeed, since contains time derivatives of the spatial metric (see (1461) and (I50p l. 
contains second time derivatives of hij. In order to find these commutators at equal time, we have to use the equation 
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of motion for h^i ,, (IB22I) . It is convenient to choose Feynman gauge and rewrite the equation of motion by subtracting 
a multiple of the trace as (c./. (IB24I) 1 


, ( 112 ) 

where is the “inverse trace-reversed” stress tensor, defined as in (IB8I) . This suggests that the nontrivial commutator 
between T^i, with cj) could potentially lead to a contribution at order k^. However, the spatial components of = 
dij,4>dy(l) -f /{D — 2), which enter , do not contain time derivatives of so there is no such contribution to 

the equal-time commutator between <!> and $; the term in (Hill) vanishes. 

Finally, we should consider a possible term in (11081) of the form 

(113) 

with similar terms appearing in (I109D . It appears that this term could contain an 0{k^) contribution. However we 
have defined $ as the value of at a point determined in terms of the geometry, so the second-order piece takes the 
form 


K^$(2)(a;) = ^V^{x)V''{x)dfj,d,,(l){x) + V^^^{x)d^(j){x), (114) 

where ^ nonlocal functional of the metric h^i, of order Then, the commutator [4>{x), h^i,{x')] is 0{k) since 

it vanishes in the absence of gravitational interactions, so the term in (111311 is 0{k^). 

We therefore conclude that (I108D and (I109D contain all the terms that can enter the commutators at order k^. 

As we have just noted, contains a time derivative of the spatial metric, and this implies a new feature of the 
algebra, as compared to the gauge-theory case: there can be a non-zero contribution to the equal-time commutators 
[$,$]. We will consider the different commutators in turn. 


1. [i>,<i>] 


In this paper, rather than working out the full structure of the commutators in detail, our main focus will be the 
question of when the commutators are nonzero, and when they become significant. We first consider commutators 
[4>, $]. These receive contributions in particular when we consider the large-mass limit, where the operators can create 
massive particles at rest. 

A simplest example of such commutators arises in the Wilson line case, where we find from (14611 . working in H = 4, 


2 poo 


[Vw^,z{x),Vw,,z{x')] = X / ds'[h^^{x + sz),hzz{x'+ s'z)] 

^ Jo Jo 

2 foo 

= -i-^<^^(x± - x'j_) / ds , 
o Jo 

using the equal-time commutator (specialized to D = 4; see (|B14p h 

[A^'"(a;),hAa(a;')] =-f 'nxa S^{x-x'), {t = t',D = 4). 


(115) 


(116) 


Thus, from (11091) . there is a nonvanishing contribution to the equal-time commutator 





dz4>{,x)dz4>{x')5'^{x 


_L — 



(117) 


where we omit terms proportional to other derivatives of (j). The divergence here can be regulated by basing the 
Wilson-line construction on a platform at a large, finite z = cutting off the integral at Z. 

We next consider commutators proportional to 4>{x)(j){x')', these are the terms that make the leading contribution 
in the large-mass, zero-momentum limit. These commutators arise from 

An interesting case is that of the commutator of two operators with the Coulomb dressing. The time derivative of 
Vq is, using (jll2l) and neglecting the stress tensor terms which we have argued do not contribute to our calculation, 
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and with r = y — x, 

v«{i) = -iL ■ '‘"'f®’’ 

47r ./ 


2r 


47r 


d^y 


VV" 1 2r 


r 


hor{y) 


47r 




Vrhrr ( 

- hauV 


2r 


2r 


+ T terms 




h i 3hrr i_i c-3/-»\ , ^Or 

- -3 -W + ^ 


<b r'^d^n 

47r 


V 7 - /Zt’ 

2r 


2r2 


+ T terms . 


(118) 


Here we have also used the identities (IC4I1 . (jCBIl . When evaluating the commutator with V(9, we can neglect the 
boundary terms. This is because the commutator of h^^{y) with Vq{x') decays as 1/r, so when we commute Vo(x') 
with the boundary term in (I118D , the resulting integrand is l/r^ and hence its integral vanishes as the surface is taken 
to infinity. 

The commutator \Vq{x),Vq{x')\ is the sum of three terms, which we denote by [ ]i, 2 . 3 . The bulk term from the 
commutator of ft-or with hor is, with = if — x', 


1 = 




327r2 


IK, 


d^ydf. 




87r|a; — x'\ 

The piece of the commutator coming from the l/r^ term in (I118|) is 

r 1 /" ,3 1 ~ 3(r • f'Y 

- ‘ d^y- 


327r2 


rSj./ 


This integral is evaluated in appendix [Cl and gives 


J2 = 


67r|a; — x'\ 


Finally, the d-function term in (I118D leads to 


Adding these terms gives 


J3 = 


967r|a; — x'\ 


(119) 


( 120 ) 


( 121 ) 


( 122 ) 


[VS{xfVSix')] = [l + [h + [h = - 


327r|a: — x'\ 


(123) 


Then, if we consider $c(a^0 to create a static source, for which <j){x') = im4>{x'), and use (11091) . this leads to a nonzero 
commutator 


Cm 

[ 4 >c(x),$c(x')] = [VS{x),VS{x'mx)f>{x') = ,^^Mx),f{x') . 


(124) 


Note the comparison between this result and the QED result (1551) : this is even clearer in the static limit (j){x) = im(f>{x). 


In higher dimensions, we expect a denominator \x — x'\^ Thus in gravity we find nonlocal commutators, and in 
particular commutators proportional to the gravitational potential. 


2. [4>, 4-] 


As noted above, a difference between gravity and QED is a nonzero contribution to the equal-time commutators 
[4>,<i)]. Such contributions arise for both the Wilson line-dressed operators, and the Coulomb-dressed operators. 
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These are most easily elucidated by considering the commutator [<l>vy^ (cc), ^(^(O)]- For simplicity, consider the case 
where x = (0, 0, z). Then, we see from (I108D that these commutators have a contribution from 


2 /•OO /» 2 

[ywA^)^^ci^)] = J dz' j z'),hij{y)]^^^ 


IK 

^2 




where in the last line we have introduced a large-z cutoff Z. Thus we find 


(125) 


[^w,(a;),d>c(0)] = / z)d^(j){x)^{Q) 


(126) 


where we don’t include terms proportional to other operators. 

While initially surprising, the log has the following explanation. The time derivative in V^, (ISU)) . creates the 
perturbation of the spatial metric which we have shown corresponds to linearized Schwarzschild, 


Khlrr - 


Wirr 


(127) 


After acting with the operator ^^(O), the proper distance from a point near infinity is thus corrected by a term 
logarithmic in r; since the Wilson-line construction uses the proper distance, this leads to the log in (11251) . Note 
that this logarithmic term is closely similar to another logarithmic dependence seen in a related physical effect: the 
Shapiro time delay. In higher dimensions, the metric perturbation created by Vq is 


Khrr oc 


Gdtti 

rD-3 ’ 


(128) 


and so the corresponding correction does not require a cutoff, and varies oc z'^~^. 

One can likewise examine the commutator ( 2 :), ‘I’Wj (a^Ol) here one finds a result that is quadratically divergent 
in the cutoff Z for H = 4. 

Finally, we consider the commutator of two of the more “physical” 4>c’s. The commutator of the dressings is 


[VS[xlV^[x')] = - 


3k^ 
327r2 
3iK^ I 


(Fy / Fy 


647r2 


Fy 


ipj 


(r • f'Y — 1/2 , 


[hjk{y),him{y')\ 




(129) 


where r = y — x and = ^ — x'; in the second line, ^ = y. The necessary integral is done in appendix [Cl yielding 


[VS{x),V^{x')] = 


iF x'* — X* 
647r lx — x'l 


(130) 


This result, then, gives a commutator 


[$c(a:), $ 0 ( 2 :')] = ^{x)di(j)[x') -f dF{x)(j>{x') 

647r L 


J X — X' 


(131) 


This term does not decay with distance, but does vanish with the momentum of the fields, and so vanishes in the 
static limit. In higher dimensions, we expect a falloff ~ |x — x'\^~^. 

Note that the commutator of the form (I131D can be eliminated by adding a linear combination of A14° and Aid* 
appearing in equations (1731) . (1741) to V/f. For example, the combination 


V^[x) = VS{x) + -/FV\x) 


(132) 


commutes with V^{x'), and so the resulting dressing does not produce a commutator (11311) . More generally, one 
expects a one parameter family of dressings with linear combinations of AV'^ and Al/* added to V/f with the same 
property. Note that these will also change the commutators but the latter can be shown to still take the 

same form as in (11231) . with different numerical coefficients. 
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C. Further comments 

We have only evaluated some indicative commutators, which reveal nonlocal behavior, as compared to LQFT, 
of our gauge-invariant, gravitationally-dressed operators. Other commutators can be likewise evaluated, with more 
effort, using similar techniques. In particular, one can evaluate the commutators of the operators with the worldline 
dressing, corresponding to more general states of motion of a particle, and find similar results to the simple cases we 
have shown. 

We also note that there are contrary claims in the literaturefl^, that gravitational dressing does not modify the 
local properties of commutators of field operators. The same claim was madefy for QED, for similar reasons. But, 
as we have detailed above, a closer inspection of the Dirac brackets for QED fixed to axial gauge explicitly shows 
noncommutativity that matches that of the dressed operators, giving one confidence in our methods and results. 
Thus, we likewise expect that the nonzero commutators for gravity - which are similar in structure to those of QED 
- are also present, and could likewise be derived through a Dirac bracket analysis. 


V. DISCUSSION 

The diffeomorphism-invariant observables that we have constructed in this paper potentially play multiple important 
roles in better understanding aspects of quantum gravity. 

A first role is to control infrared divergences in scattering. While that has not been a primary focus of this paper, 
we have noted that the worldline-dressed operators described in section III Cl have been argued by [sl, HH to regulate 
IR divergences in scattering in QED.^ Section [III Cl has constructed analogous operators in gravity. Thus, we expect 
the corresponding gravitational dressings to analogously address IR divergences in gravitational scattering. We leave 
development of such a treatment of scattering to future work; for related work see [i^ . 

Another very important role is in capturing features of the fundamental structure of quantum gravity. As we 
have noted, locality in LQFT is most clearly described in terms of commutativity of subalgebras of gauge-invariant 
observables associated with spacelike-separated regions of spacetime. A key question for quantum gravity, if it is a 
quantum-mechanical theory, is thus what algebraic structure governs its observables, and how this structure relates 
to possible localization, and reduces to the locality structure of LQFT in the limit G —>■ o[^. 

While the strong gravitational regime still poses many puzzles, we have found that even the weak-field regime allows 
us to infer apparently important aspects of this algebraic structure, since non-trivial results can be found perturbatively 
in G, and since any more complete algebra relevant to arbitrarily strong fields must have these contributions at leading- 
order in an expansion in G. In particular, we have found, at order G in an expansion of commutators of the gauge- 
invariant observables, obstructions to commutativity of operators associated to different “regions” of spacetime. If 
locality is defined in terms of such commutativity, it fails for the gravitationally-dressed operators we have considered, 
and in a way that suggests that it will not necessarily be easy to restore with definitions of different operators. This 
appears to be an important structural aspect of a quantum theory of gravity, if it is to agree with quantized general 
relativity in the weak-field regime. 

One key piece of information is that of when the noncommutativity becomes significant; that is expected to be 
an important characteristic of the “correspondence boundary” where quantum gravity reduces to LQFT [11113. For 
example, for two particles of energy A, (11241) indicates that the dressing only makes a small correction to <I>(x)<I>(a;') 
when GE < \x — x'\, or, for general D, for 


GdE <\x-x'f~^ . (133) 

This is in accord with the locality bound proposed in which stated that LQFT ceases to give an accurate 

description of the state once the bound (I133D is violated. Note that we have also found significant corrections in 
the $(a;')] commutator, which appear to become relevant in a regime where GpE ^ — x'\^~‘^, where p 

is a characteristic momentum of the particles being created. While interesting, and clearly related to the physical 
effect described in (11261) , these commutators can be removed by modifying the dressing, as described in the preceding 
section. 

Notice that we can, at the order G to which we work, achieve commutativity of the gravitational Wilson line 
operators of (HSl) associated to spacelike separated points x and x', if we run their associate Wilson lines in different 
directions so that they also stay spacelike separated. However, such operators do not have a clear identification with 


Connections to and issues with the related analysis of |41| are discussed in 1421 . 
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a compact spacetime region, and appear to be more clearly associated with a noncompact neighborhood, extending 
to infinity and containing the flux line. Moreover, we expect important correctionsat higher order in G. First, the 
infinitesimally-thin Wilson lines have infinite energy density, and thus are expected to receive significant corrections 
once self-coupling of gravity is taken into account. This is expected to “thicken” these Wilson lines. Moreover, it 
would appear that the thickness of the corresponding region would grow as the mass or energy sourcing the flux 
lines increases. This suggests an important lesson for any putative subalgebra structure: a monomial in the operator 
(^) would appear to be associated with a larger and larger region as the order of the monomial grows; moreover, 
this also strongly suggests that such a monomial doesn’t commute with for high enough order, spoiling any 

commutative subalgebra structure. 

While we have treated the case of asymptotically Minkowski space, many of the features we have described should 
carry over to the case of anti de Sitter space with relatively minor modification. In particular, one can likewise construct 
Wilson line operators, associated with “Fefferman-Graham” gauge, there; for related discussion see [ll,[ll,[3^. These 
can alternately be averaged over directions, as we have done in section IIIIBl to produce a Coulomb-like dressing. 
One would then find similar commutator structure for the corresponding operators. In particular, two Wilson line 
operators with overlapping Wilson lines are not expected to commute, as in (I93p . This differs from a claim of [l8l| . 
though we have understood the origin of the conflict in the simple example of QED in section IIV A 31 Also note that 
for many purposes working in AdS will provide an effective infrared regulator, with characteristic length ~ i?AdSj to 
calculations with the flat-space operators described in this paper. 

The present discussion also has interesting relations with the similar problem of observables in de Sitter space - 
though here nontrivial features are encountered. In particular, [1, [l^ discussed formulation of observables where field 
operators are separated by a given geodesic distance; as in the present work, these will create particles together with 
a gravitational flux line connecting them. Thus, for far separated particles, or if one particle is taken to be massive 
and provide a “platform” with respect to which we measure, the construction for the remaining particle is very similar 
to the ones we have described. However, notice that in such a picture field lines appear to terminate on the pair 
of particles; none can reach asymptotic infinity since space is compact. This new feature arises from the nontrivial 
nature of the de Sitter background. 

Past examples^ of diffeomorphism-invariant, approximately local observables such as those just mentioned are 
relational, in that the position at which we create or measure a particle is defined in relation to other particles or 
features of the background state. Interestingly, the observables constructed in this paper do not require any such 
local structure to define location; the location of the field operator is defined in relation to structure at infinity. Thus 
these are still relational, though in a somewhat different fashion. 

A final role to consider for such observables is their connection to observation or experiment. The observables of this 
paper, like those of @,[11 are observables in the usual mathematical sense of quantum mechanics - they are Hermitian, 
gauge-invariant operators on the Hilbert space. However, they do not have an a priori connection to observations 
made by “observers inside the system” and thus were referred to as q-observables in [dl. Some such q-observables 
are expected to be related to observations such observers can make (thus to what experimental physicists would call 
“observables”); further development of this story is left for future work. 
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Appendix A: QED basics 

Here we collect some basic formulas relevant to quantization of QED. 

The Lagrangian of QED takes the form 

Cqeu = + Fm ■ (Al) 

The second term is a “gauge-fixing” (really, gauge-invariance breaking) term; gauge transformations act as 

A^(x) A^(x) - d^A(x) . (A2) 

Gauge symmetry is restored for a = oo, a = 1 gives “Eeynman gauge,” and a —>■ 0 gives “Lorenz” or “Landau gauge.” 
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The third term is the matter Lagrangian. The corresponding equations of motion are 


1 




a 


1 


where f d^xCm is the current; here □ = d^d^. 

A particular matter Lagrangian is that for a charged scalar, 

with gauge transformation 


a 


Cm = -\DM^ - 


(t>{x) 


^—iqA.{x 


Uix). 


and with covariant derivative 


(A3) 


(A4) 

(A5) 


The corresponding current is 


The canonical momenta are 


D^4> = - iqA^fj) . 


(A6) 

(A7) 


T,^ = {Do4'r , = (AS) 

and 

TT* = = -E\ 7r° = -duAf^ . (A9) 

a 

At a = oo, 7r° of course vanishes, yielding a constraint. The equal-time commutators are 

['K^{x),(j}ix')]t =-i6^~\x-x') = [TTl{x),(l)*{x')]t (AlO) 

\t=t' \t=t' 

and 


[k^{x), Ay{x')\\ =—id^S^ ^{x — x'). (All) 

I t = t' 

These commutation relations provide initial data for the unequal-time commutators. In the free limit, these satisfy 
the free equations of motion. For the free scalar, we have 

[(l){x),(j)* {x')] = iA{x — x'). (A12) 

where A is the massive Pauli-Jordan function, satisfying 

(□ — m^)A(x) = 0 , A(x)| =0, 9tA(x)| =—S^{x). (A13) 

Note also that 


A(x — x') = Ga(x, x) — Gr{x, x) , (Al4) 

where Gr and Ga are the retarded and advanced Green functions, respectively, 

Gr{x,x') = i9{t - t'){Q\[(j){x),(j)*{x')\\Q) , Ga{x,x') =- t){Q\[(j){x),(j)*{x')\\Q) . (A15) 

The commutators for the electromagnetic field may be written in terms of the massless Pauli-Jordan function (again, 
advanced minus retarded Green function), satisfying 

UD{x) = 0 , (A16) 

D{-x) = -D{x) , (A17) 

£>(x)| =0, (A18) 

1 1=0 

dtD{x)\ =—5^~^{x). (A19) 

11=0 


25 


In 4 dimensions, D{x) is given by 


where e{t) is the sign function. The commutators become 

[A^{x),A^{x')] = iD^^{x - x') 


(A20) 


(A21) 


where 


D^,u{x) = VtJ.i^D{x) + (1 - a)d^d^E{x), (A22) 

and OE = —D. More explicitly, E{x) is the Green function for the operator with the boundary conditions 

E{x)t =dtE{x)< =dfE{x)< =0, dfE{x)< =-S^-\x) . (A23) 

11=0 11=0 I t = 0 11=0 

In D = 4, it is constant in the forward lightcone, and in the backward lightcone, and is given explicitly by: 

E{x) = -^em-x^), (A24) 

Eq. (IA22|) can be verified by checking that the resulting commutator satisfies the equation of motion, 

+ (^ - l) 

with initial conditions given by: 

[AQ{x),AQ{x)\ = —ia5^~"^(x — x) , [Aiix), Ajix')] = i5ij5^~^{x — x) , (A26) 

which follow from (lAIIl) . 

Note that the quantity B = dfj,A^ generates gauge transformations on both the electromagnetic and matter fields: 

i[B{x'), Af^{x)] = adfj_D{x'— x) , i[B{x'), (j){x)] = iaqD{x'— x)(j){x). (A27) 

This is an infinitesimal gauge transformation, 

i[B{x'),Afj^{x)] = dAAf,{x) =-dfj^A{x) , i[B{x'),(f>{x)] = Sa^x) =-iqA{x)(j){x') (A28) 

with A{x) = aD{x — x'). Here the commutators with follow from (IA2II) . To check the 4> commutator, we first use 

the equal-time commutation relations to show that it holds at equal times, 

[B{x'),(l){x)]\ =0 , [d[jB{x'),(j){x)\\ = —aqS^ix — x')(j){x) , (A29) 

\t=t' \t=t' 

where we have used the equation of motion to write d[jB = a{diTr'^ + The identity (IA27|) then follows at unequal 
times using the fact that B satisfies the free equation of motion DH = 0 when the current j'^ is conserved. 

Thus, a gauge-invariant operator constructed out of the electromagnetic field and the scalar field (j) will commute 
with B{x) . The physical states of the theory are those annihilated by the positive-frequency part of B: 

B+{x)\ijj)=0 \^p) is a physical state. (A30) 

Then if $ is a gauge-invariant operator, it will commute with the positive-frequncy part of H, and hence maps physical 
states to physical states. 


D^„{x) = 0 , 


(A25) 


Appendix B: Gravity basics 

Here we collect some basic formulas relevant to perturbative quantization of gravity. 

The scalar Lagrangian density for Einstein gravity takes the form 

•^grav = + Egf + Cm , (Bl) 
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where = 327rG'D, Gd is the D-dimensional Newton’s constant, R is the scalar curvature and Cm is the matter 
Lagrangian; a particular example is that for a scalar with mass m, 

= • (B2) 

The second term in (jBlI) . Cgf, denotes a “gauge-fixing” (really, gauge-invariance breaking) term. If one picks a 
background metric g^, one useful choice is 


VWl^gf = - 


1 

|g| 




(B3) 


where denotes the covariant derivative with respect to g^ and g^ is used for the contraction of the v index. Gauge 
symmetry is restored for a = oo, a = 1 is “Feynman gauge,” and a 0 is an analog of “Landau gauge,” which 
enforces the de Bonder gauge condition 


v° = 0 . 


(B4) 


When g^ is the flat metric, this reduces to the usual harmonic gauge condition, which can be expressed in any of the 
equivalent forms: 


{VWlg^'') = 0 , = 0 , = 0 . (B5) 

The latter says that the coordinates are harmonic functions of spacetime. 

For the purposes of this paper we primarily focus on the linearization of gravity about flat space. The metric 
perturbation is defined by 


. (B6) 

For the linearized theory we need the quadratic-order expansion of ^/\^\Ceh = 2\/[^i?/K^ in h. Various formulas 
simplify if we define the “trace-reversed” metric perturbation, 

hfiu — ; (B*^) 

with h = the inverse to trace reversal, in D spacetime dimensions, is 

^/ii/ — ^/ii/ = h^y . (B8) 

The quadratic part of the Einstein-Hilbert action can then be simplified to 

= -^d^h^yd^h>^‘' +d^hx^dyh''>^ +t.d. (B9) 

where indices are raised with the fiat metric 77 and the last term is a total derivative. Likewise, to quadratic order, 
the gauge-fixing term gives 


1 

' 2 a 

and so the combined quadratic action for gravity is (dropping total derivatives) 




{CeH + Cgf) = -id^h^yd^h^'' + ( 1 - - ) ■ 


1 


(BIO) 


(Bll) 


This action exhibits the simplicity of Feynman gauge, a = 1. Here, the canonical conjugate to hg^y is so we 
have the equal-time commutation relations 


[hgy{x),h^’^[x')]\ =i5^16^ ^{x-x'), 


(B12) 
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where we define 


5^Z = 


(B13) 


with symmetrization convention + A'^^)j2. Equivalently without using the trace-reversed field, 


I t=t' 



D-2 J 




f'). 


(B14) 


When K = 0, the field equation reduces to = 0, from which we find the unequal-time Feynman-gauge 

commutation relations for fields in the interaction picture: 

[hf,^{x), h^‘^{x')] = i5lZD{x - x') , (B15) 


or 


[h^,{x),h>'^{x')] = I D{x - x') , 

with D{x) given in appendix!^ Alternately, the momentum-space two-point function is 

{p)) = -^S^^{2Tr)^S^{p + p') , 

with the corresponding correlators for the metric found via the transformation (IB8I) . 
These expressions can be generalized to a ^ 1; let us introduce the variable 

^=1--. 
a 


Then, the quadratic action (jBllIl takes the form 


\/\F\ {CeH + Cgf) 




Act 


where Lg^^'^{l3) is the second-order linear operator defined by 


(B16) 


(B17) 


(B18) 


(B19) 


A(7 


-'flu 



1 


X<T 


d^dp 


2/3 






Defining the stress tensor as 


(B20) 




(B21) 


the linearized gravitational equations then take the form 

Lg,^^/3)hxa = , (B22) 

where in this equation the metric in Tpi, gets replaced with rj to leading order. Then, the propagator for the metric 
takes the form 


{Thp^{x)h^'^{x')) =i{L ^)p^^'^{x,x') . 
Note that the linearized gravitational equations (jB22l) simplify to 

ni -_ 1 t 

'—'^flU flU 


(B23) 


(B24) 


in Feynman gauge (/3 = 0). 
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One can also work out the canonical momenta and commutators for the metric. From (IBllI) . one finds 

Tr°i^ = fi^i^ - , 

_ i3d^h'^°S^^ . (B25) 


When /3 = 1 there is no gauge fixing, so the momenta do not contain time derivatives, and give the expected 
constraints. For general j3 the canonical commutators are 

['K>^'^{x),hxa{x')]\ =-i5^l5°~^{x - x') . (B26) 

I t = t' 

When /? = 0, this agrees with the Feynman gauge result (IB12I) . 

These commutation relations provide initial data for the unequal-time commutators, which in the free limit satisfy 
the free equations of motion, (IB22I) with T = 0. These are 

[h^yx), h^'^{x')] = iD^^^'^{x - x') (B27) 


with 


D^,^<^(x} = S;::D(x) + 2(l-a) 




(m 


E{x) 


(B28) 


and D{x) and E{x) given in appendix [Al This can be checked by verifying that 

Ll^J^'^Dxa^\x,x')=0 (B29) 

and that satisfies the initial conditions implied by 

[h^^{x),hxa{x')]\ =0 (B30) 

I t = t' 

and (|B26I) . 

The gauge condition b'^ = generates infinitesimal diffeomorphisms, just as the gauge condition in QED 

generates infinitesimal gauge transformations (c./. (IA28I) 1. From (IB27I) and (IB28I) we find the commutator 


i[b'^{x'),h^yx)] = -d^yy\x) - d„yy\x) , (B31) 

with the infinitesimal diffeomorphism generated by the vector fields 

^I^H^) = ^S<;,D{x-x'). (B32) 

One can also see that 

i[b'^{x'),(l){x)] = -Kyy\x)df^(l){x) . (B33) 

To show this, one first checks the equal-time versions of it. 


i[b'^{x'),(j)ix)]t = 0 , 

I t=t' 

ir{x')A{x)]\ ^=i^[T°‘^{x')A{x)] = -^d^Hx)S^-yx-x') (B34) 

where the latter commutator follows from the equations of motion (IB22I) . Then, by taking the divergence of the same 
equations of motion, we find 


ab‘^ = 0 (B35) 

when the stress tensor is conserved. Thus since (IB33|) satisfies this equation of motion in x' and the initial conditions 
(IB34I) at t = t', it holds for all x'. 

In particular, this implies that if $ is any operator invariant under linearized diffeomorphisms, then the linearized 
metric h defined to leading order in k by 

h^i„{x)<^{x') = [h^yx),^{x')] (B36) 

satisfies the gauge condition = 0. 























29 


Appendix C: Some useful formulas 


Here we collect some derivative formulas and integrals used in the main text. It will be useful to introduce a radial 
vector field such that = 9^.; we will denote the spatial metric by qap Then we have the identities (given for 
H = 4, though easily generalized to D 4) 




= - 
r 


rafi3\ _ q^iaVp + qtipTa “ 3 f^faf ;3 


□ 


f^fp 

r 

fgfp 


rp 


= 2 


qa/3 - ‘jfgrp 47r 




dadpd^r = —{fafpf^ - r(aqp-t)) 

dadpdj- = - 4 - ^q{apd^)S [x) , 

T T 0 

where in the last two equations fi^aqp-y) = {faqp'y + rpq~fa + f-yqap)li- 

Calculating the commutators in section ITVl requires evaluation of certain integrals. The first is 


(Cl) 

(C2) 

(C3) 

(C4) 

(C5) 

(C 6 ) 

(C7) 

(C 8 ) 


(C9) 


where r = y — x, f = y — x'. This is evaluated by choosing spherical coordinates based at a; = 0, and with polar 
direction defined as that of d = a;' — a;. In particular, we then find 


(f • f') =1-^(1 —cos"^0). 


(CIO) 


We then expand 1 /r' and 1 in Legendre and Gegenbauer polynomials respectively: 

1 1 [^Sz^oa(cosd)(^)^ d<r 

x' \/r 2 — 2dr cos 9 + cP \ 3 ^(cosd) r<d 

1 _ _ 1 __ / 7 ^E“oC'f^^(cos 6 l) (^)' d<r 


j ,/3 (r^ — 2 drcos 0 + 


E“oC'f''^(cosd) ( 3 )' r<d 


The angular integral picks out the I = 0 term of each sum, and we find 


d^y 


1 — 3(f • f') 


>\2 


= 27r y r^drdcos(d) 


-2 


3d' 


7 + 33^(l-cosd ) 


= 27r 


dr — 


4 4 \ r 

rd rd) X 


dr ( - — 


4d2 

TiT 


IOtt 


3|a; — x'\ 


(Cll) 

(C12) 


(C13) 


Another needed integral is 


r = 


3 a(f • f'y + b 
d-^y ^ -r'* 


(CM) 




















30 


where a, b are fixed constants. For this we center spherical coordinates about the point x', with the polar axis 
determined by d = x — x'. By symmetry, the integral must be proportional to d\ with coefficient d^P. We also use 
(|C10p . with giving 


/ dr'd cos 9 

a + b d^(l —cos^0) 

a 3 

J 



cos 9 . 


(C15) 


As before, we expand 1/r and 1/r^ in Legendre and Gegenbauer polynomials, using (ICllI) . (IC12I) with r r'. In 
both cases the 9 integral picks out the I = 1 term of the sum: 

2 4 

dcos9 cos9Pi(cos9) = -dll, / dcos9cos9C^^^(cos9)(l — cos9^) =-dn. (C16) 

-1 3 5 


Substituting this into the above integral, we find 

nd 


ddP = 27r 


/o 


o ~ ^ o + ^)i72 “ 7® 


and thus 


' 3 a(f ■ f'Y -n o 
d?y- --r'* = 27 r 


4 d? 


5 r" 


= 27r 






(C17) 
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